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Multidimensional persistence behaviour in an Ising system 
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We consider a periodic Ising chain with nearest-neighbour and r-th neighbour interaction and 
quench it from infinite temperature to zero temperature. The persistence probability P(t), measured 
as the probability that a spin remains unflipped upto time t, is studied by computer simulation for 
suitable values of r. We observe that as time progresses, P(t) first decays as t~ ' 22 (-the first regime) , 
then the Pit) — t curve has a small slope (in log-log scale) for some time (-the second regime) and at 
last it decays nearly as t -3//8 (-the third regime). We argue that in the first regime, the persistence 
behaviour is the usual one for a two-dimensional system, in the second regime it is like that of a 
non-interacting ('zero-dimensional') system and in the third regime the persistence behaviour is like 
that of a one dimensional Ising model. We also provide explanations for such behaviour. 

PACS numbers: 64.60.Ht, 05.50.+q 



I. INTRODUCTION 



The tendency of a spin in a spin-i Ising system to re- 
main in its original state following a quench from infinite 
temperature to zero temperature has been extensively 
studied over the last decade and is an example of the 
phenomenon called persistence in dynamical systems [l|- 
@. The probability P(t) that a spin does not flip upto 
time t, exhibits a power law behaviour 



P{t) ~ t~ 



(1) 



where 9 is a non-trivial exponent, as it is not related to 
any other static or dynamic exponent. In one dimension 
with nearest-neighbour interaction, it has been proved 
exactly that this exponent is = 3/8 |3j. In two di- 
mensions, again with nearest-neighbour interaction, this 
exponent has been numerically evaluated [6( as — 0.209. 
For Glauber dynamics, the one-dimensional Ising system 
is equivalent to a one-dimensional A + A — > diffusion 
system (see below). Hence, the "zero-dimensional" per- 
sistent behaviour should correspond to that of a system of 
non-interacting particles initially spread randomly over a 
chain with density p and then diffusing independent of 
each other. (The rule of diffusion is to take a step to the 
right or to the left with probability 1/2.) The persistence 
probability for this system has been shown Q to decay 
stretched-exponentially, 




Pit) = exp[ 



(2) 



The objective of this communication is to report the ob- 
servation that an Ising system with nearest-neighbour in- 
teraction on a rectangular helical lattice, when quenched 
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FIG. 1: 

The lattice with helical boundary condition. 



to zero-temperature from an infinite temperature shows 
zero-, one- and two-dimensional persistence behaviour in 
different regions of its temporal evolution. In the next 
section, we shall describe the details of the system simu- 
lated and the algorithm followed, alongwith the results. 
An explanation of the simulation results will also be pre- 
sented in Sec. III. In the last section we shall discuss 
some subtle issues. 



II. THE ALGORITHM AND THE RESULTS 

In this section we shall first describe the details of our 
system and our algorithm, and then present the results. 
Consider Ising spins s\, S2, ■ ■ ■ , sl placed on a chain of 
L sites with the periodic boundary condition i + L = i. 
Initially, each spin is +1 or —1 with equal probability. 
One iteration of the system consists of the following 4 
steps: 

(i) choose one (say, fc-th) spin randomly, 
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FIG. 2: 

Plot of P(t) as a function of t for r = 20. The curve is for 
L = 10000 and averaged over 100 configurations. The first, 

second and third regime are marked in the figure. The 
curves for L — 5000 and L = 20000 fall on the curve shown 

here. 



(ii) calculate its energy 

E k = Sfc(Sfc_i + Sfc+l + Sfc-r + S k +r), (3) 

(iii) flip Sfc with probability 1 if E k > 0, and with prob- 
ability \ if E k = (do not flip at all if E k < 0), 

(iv) repeat the steps (i) to (iii) (L — 1) times more (ran- 
dom updating). 

Here, r is a parameter of the model and must lie between 
2 and L. Clearly, we have nearest and r-th neighbour in- 
teraction and our system is equivalent to a nearest neigh- 
bour rectangular Ising system of size r x (L/r). The lat- 
tice is not periodic in the two axial directions, rather the 
chain is wound as a helix of periodicity r with the ends 
(first and L-th sites) put side by side (Fig. [T]). When 
r <C VL, the system is effectively a one-dimensional one, 
while for r ~ yL it is a two-dimensional one with as- 
pect ratio a — r 2 /L. In this communication we study 
the case of a ~ 0.01, so that the system is effectively 
one-dimensional, or at least a narrow strip. We do not 
consider the situation when the condition r -C \TL is not 
satisfied. 

After every iteration we compute the fraction of spins 
that has not been flipped till now. This fraction, aver- 
aged over many realisations of the system gives us the 
probability of persistence P{t). We present in Figs. [2] 
and [3] the simulation results for this quantity. At first 
Pit) decreases as t~ e with 9 = 0.22 ± 0.01 showing two- 
dimensional behaviour. This behaviour continues upto a 
certain time, say, t%. For a given L, the value of T\ has 
been observed to increase with r as r 2 (Fig. 2]). We call 
the region 1 < t < t\ the first regime. Next follows the 
second regime extending upto some T2 iteration where the 
Pit) — t curve has a small slope in the log-log scale. We 
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FIG. 3: 

Plot of P(t) as a function of t for r = 10, 20 and 30. The 
curve is for L = 10000 and averaged over 100 configurations. 
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FIG. 4: 

Variation of ri vs r. The curve is for L = 10000 and 
averaged over 100 configurations. 



shall see below that in this regime the system behaves as 
a zero-dimensional one. For limitation of computational 
resources, we could not achieve precise evaluation of ti, 
but could observe that, like ti, this quantity increases 
with r. At last comes the third regime (for r 2 < t < oo) 
where P(t) decays nearly as t with 9 = 0.375 ± 0.01 
showing the one-dimensional behaviour. 

III. EXPLANATION OF THE RESULTS 

We now explain the observations in the three regimes 
one by one, by approximate analytic arguments but an 
exact analytic calculation of the persistence behaviour 
spanning over the three regimes is yet to be done. 

The First Regime : Here the system shows normal two- 
dimensional behaviour. The persistence curve P{t) vs. t 
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FIG. 5: 

Presence of only a small number of kinks in the domain wall 
(A particles) in the second regime. Dynamics occurs only at 
the kinks. 



saturates at a time t\ ~ r 2 to a value that varies as r~ zB 
where z is the dynamical critical exponent (~ 2) [J S] . 
After the system reaches the saturation stage, the domain 
walls (lines separating unlike spins) are mostly parallel to 
the X axis, with only a small number of kinks (Fig. [5j. 
The dynamics occurs only at these kinks, which we call 
l A particles'. For the updating rule stated above, each 
A particle jumps to the left or right with probabilities 
1/2, 1/2, assuming that the walls are sparsely distributed 
over the system (an A particle does not have another at 
a distance < r) . The density pa of A particles, measured 
as the number of vertical domain walls per site is shown 
in Fig. [S] This quantity also shows a plateau region in 
the second regime, like persistence. It is important to 
note that, for a fully periodic system shown in Fig. [3 
there can be only an even number of A particles in a 
row. Presence of one particle in a row is hence ruled out 
and two or more particles get annihilated within time n . 
The dynamics therefore stops completely at t = t\ in the 
case of a fully periodic system, and the second and the 
third regimes do not appear. 

The Second Regime : What happens to our (helical) 
system for t > t\ ? For an interval of time t\ < t < T2, 
the A particles (kinks in the domain walls) are so far 
away from each other that they cannot "see" each other 
and diffuse independently. For a system of independent 
random walkers of density p, the persistence has been 
shown 0] to obey Eq. To compare our persistence data 
with Eq. we assign each site to be persistent at t = t\ . 
This makes P{t\) — 1 and obliterates the distribution 
of persistent sites created in the first regime. Then we 
note down the (almost constant) density pa in the second 
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FIG. 6: 

Plot of density pA of A particles as a function of t for r = 10 
and 20, L = 10000. The data was averaged over 1000 
configurations. In the nearly horizontal region (the second 
regime) p A = 5.50 x 10~ 3 for r = 10 and p A = 1.44 x 10" 3 
for r = 20. 




FIG. 7: 

The lattice with periodic boundary condition. 



regime and observe that the P(t) data here obeys the 
relation 



P{t) = 1 - ap A Vt^Ti 



(4) 



with a w 1.6 (Fig. [8]). Since the value of the slope a 
is close to 2y/2/T: — 1.596, and since the value of pa is 
small, Eq.[5]is obeyed and we conclude that in the second 
regime, the system behaves as one of zero-dimension. 

The Third Regime : When the system evolves further 
(t > r 2 ), the A particles start meeting (and annihilating) 
each other and the usual one-dimensional dynamics leads 
to t~ e behaviour with 9 « 3/8. Since the density pa 
is too small, one needs to average over a large number 
of realisations and simulation of the third regime is a 



4 



D- 

0.94 - 
0.92 - 




0.9 1 1 1 1 1 1 1 1 1 1 1 

5 10 15 20 25 30 35 40 45 50 

(t-x,) " 5 
FIG. 8: 

Plot of P(t) as a function of y/t — tl in the second regime of 
the persistence curve for r = 20 and L = 10000. P(t) was set 
to be 1 at ti — 500. The curve fits to 1 — 1.58pAy/t ~ T i 
with pa = 1-44 x 10 -3 . The data were averaged over 1000 
configurations. 



computationally intensive job. The P(t)—t curves remain 
the same for different values of L for a given r, but get 
shifted (maintaining £~ 3 / 8 behaviour) as one varies r at a 
given L (Figs. [Hand [3]). That the dynamics in the second 
and the third regimes is indeed described by simple one- 
dimensional A + A — * dynamics is further corroborated 
by two numerical experiments : 

(i) If we turn off the r-th neighbour interaction at t = t\ , 
the slope of the persistence curve (in log-log scale) does 
not change much (Fig. [9]), indicating that it is chiefly the 
nearest-neighbour interaction that drives the dynamics. 

(ii) Let us consider a periodic chain of L sites and sprinkle 
randomly some particles (excluding multiple occupancy 
at a site) with density p. Starting with a low (~ 0.005) 
value of p, we let the system evolve according to the 
usual A + A — > dynamics. The result is the (zero- 
dimensional) second regime (Fig. [TO]) for the first 100 
iterations, followed by the (one-dimensional) third regime 
(Fig. QT]) . The second regime is found to follow Eq. [2]with 
a = 1.55 and the third regime shows the usual i~ 3 / 8 
behaviour. 



IV. DISCUSSION 

(1) We shall first discuss a subtle issue regarding the 
dynamics in the third regime of the Ising model. When 
two "kinks" of Fig. [5] come closer than r to each other, 
a domain of length less than r is formed. The rules of 
update renders (i) every spin lying within this domain 
liable to flip (with probability 1/2) and (ii) the two spins 
at the ends of this domain bound to flip (with probabil- 
ity 1). The dynamics thus differs from the one for the 
A + A — > model, since once the domain is less than 
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FIG. 9: 

Plot of P(t) as a function of t for simultaneous nearest 
neighbor interaction and long-range interaction r = 20. 
After 100 iterations the long-range interaction is switched off 
(continuous line). The curve is for L = 10000 and averaged 
over 1000 configurations. The dotted line corresponds to the 
usual case, when the r-th neighbour interaction continues for 
the entire range of time. 
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FIG. 10: 

Plot of P(t) as a function of \ft for A + A — > dynamics in 
log-linear scale. At t = 0, the A particles were randomly 
spaced with density pA = 0.005. Here L = 4000 and the 
results were averaged over 200 realisations. The initial 
portion (t < 400) fits to Eq.EJwith a = 1.55. 



r in length, it is swapped within the next r steps. But 
since we are primarily interested in the region f>r, this 
difference is not of much consequence. One can flip the 
spins only at the ends of the domain (keeping the ones 
within the domain unflipped) if we replace the energy 
expression of Eq. [3] by 

E k = S k [s k -1 + Sfe+l + K(sk-r + Sk+r)} (5) 

and choose k < 1. We have checked that the persistence 
behaviour still remains almost the same. 

(2) In the case of A + A — > model with initial den- 
sity pa = 0.005, we have shown a i -3 / 8 behaviour in 
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FIG. 11: 

Plot of P(t) as a function of t for A + A — *• dynamics in 
log-log scale. At t — 0, the A particles were randomly spaced 
with density pa = 0.005. Here L — 4000 and the results 
were averaged over 100 realisations. The algebraic portion 
fits to 13.4£-°- 34 



Fig-EJ However, one comes across anomalous behaviour 
at somewhat larger system size (Fig.[T2]). Thus, after the 



third regime continues for some time, one comes across 
another saturation region, with the saturation value in- 
creasing with increase of L. We did not observe such 
anomalous behaviour for the Ising system. Work is in 
progress on this issue. 

(3) It is interesting to note that when the range r of 
the long-range interaction varies from site to site and 
takes any integer value chosen randomly between 2 and 
L, the persistence behaviour almost vanishes 9]. Since 
we have found the persistence behaviour to remain upto 
time Ti ~ r 2 , for a given value of r we conclude that 
the random nature of r tends to remove the persistence 
behaviour. 
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FIG. 12: 

Plot of P(t) as a function of t for A + A — > dynamics in 
log-log scale. At t = 0, the ^4 particles were randomly spaced 

with density pA = 0.005. The size of the samples are 
indicated in the figure and the results were averaged over 50 

to 200 realisations. The saturation region rises with 
increasing system size. Just before saturation, the curve does 
not show algebraic decay over any appreciably long region. 




